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Scanning tunneling microscopy images of hydrogen-terminated Si(110) surfaces are studied using 
first-principles calculations. Our results show that the calculated filled-state images and local den- 
sity of states are consistent with recent experimental results, and the empty-state images appear 
significantly different from the filled-state ones. To elucidate the origin of this difference, we ex- 
amined in detail the local density of states, which affects the images, and found that the bonding 
and antibonding states of surface silicon atoms largely affect the difference between the filled- and 
empty-state images. 

PACS numbers: 73.20.At, 68.37.Ef, 71.15.Mb, 73.20.-r 
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As the miniaturization and high integration of elec- 
tronic devices progress, the atomic and electronic struc- 
tures of silicon surfaces become areas of intensive re- 
search. Scanning tunneling microscopy (STM) [J has 
provided much useful information on surface configura- 
tions with an atomic resolution. Many studies revealed 
that surprisingly complex structures occur in seemingly 
simple processes such as the adsorption of atoms on semi- 
conductor surfaces. First-principles calculations are also 
implemented to characterize surface atomic configura- 
tions and to interpret electronic structures observed by 
STM. So far, there have been many studies that reveal 
interesting atomic and electronic structures on semicon- 
ductor surfaces by the combination of experimental and 
theoretical procedures 0. 

Nowadays, a Si(110) surface has been considered to 
be of basic importance in surface science and semicon- 
ductor technologies, because the hole mobility of Si(110) 
with radical oxidation is achieved to be 2.4 times higher 
than that of Si(001) 0. However, the atomic and elec- 
tronic structures of Si(110) surface have not been studied 
intensively because it is difficult to obtain atomically flat 
Si(110) surface in experiments. Recently, Arima et al. 
observed atomically resolved STM images of hydrogen- 
terminated Si(110) surfaces after wet cleaning and found 
certain characteristic features such as a zigzag chain in- 
side a terrace, a single-row bright spot at a step edge and 
an isolated zigzag chain on a terrace when the sample is 
rinsed with ultrapure water 
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In this study, we carry out first-principles calculations 
based on the density-functional theory 0, |(| to investi- 
gate electronic structures of various hydrogen-terminated 
Si(llO) surfaces. Our results show that the calculated 
STM images of filled states and local density of states 
(LDOS) are consistent with the results obtained by ex- 
periments 4,7]. On the contrary, the empty-state images 
markedly change from the filled-state ones. The zigzag 
chains run on the silicon-atom rows of the topmost layers, 
while they are formed between the silicon-atom rows in 
the filled-state images. Moreover, the width of the trough 
composed of a missing monohydride row is greater than 
that in the filled-state images, and the contribution of the 
isolated monohydride row to the images is small. These 
results are explained using the charge distribution of the 
bonding and antibonding states involved in tunneling. 

Our first-principles simulation is based on the real- 
space finite-difference method @,IUE1E|]> which enables 
us to determine self-consistent electronic ground states 
and optimized atomic geometries with a high degree 
of accuracy, using the timesaving double-grid technique 
|ToL ITU Il2l| an d the direct minimization of the energy 
functional [13|. The norm-conserving pscudopotentials 
[l4j of Troullier and Martins are employed to describe 
the electron-ion interaction. Exchange-correlation effects 
are treated with the local density approximation [l6| 
in the density- functional theory [a, |6j- The nine-point 
finite-difference formula (i.e., the formula with N=4) |8| is 
adopted for the derivative arising from the kinetic-energy 
operator. We take a cutoff energy of 23 Ry, which corre- 
sponds to a grid spacing of 0.65 a.u., and a higher cutoff 
energy of 210 Ry in the vicinity of the nuclei with the 
augmentation of double-grid points jToL ITTl ll2T | . 
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FIG. 1: (Color) Top views of topmost and second layers of 
H-terminated Si(llO) surfaces: (a) terrace model, (b) step 
model, (c) isolated monohydride row model, and (d) missing 
monohydride row model. Pink and blue balls represent H and 
Si atoms, respectively, and atoms are denoted by large and 
small balls according to the distance from the surface. 



Figure ^ depicts the top views of the computational 
models adopted here. We employ a technique that in- 
volves the use of a supercell whose size is chosen to be 
L a =61.58 a.u., L y —7.27 a.u., and L 2 =40.00 a.u., where 
L x , L y , and L z are the lengths of the supercell in the 
x, y, and z directions, respectively. Here, the direction 
perpendicular to the surface is chosen to be the z di- 
rection. To eliminate completely unfavorable effects of 
atoms in neighboring cells that are artificially repeated 
in the case of the periodic boundary condition, the non- 
periodic boundary condition of vanishing wave functions 
out of the supercell is imposed on the z direction, while 
the periodic boundary condition is adopted in the x and 
y directions. The supercell contains five silicon layers 
(i.e., the thin- film model). The topmost layers contain 
12, 6, 2, and 10 Si atoms for the (a) terrace, (b) step, (c) 
isolated monohydride row, and (d) missing monohydride 
row models in Fig. ^ respectively. The silicon atoms 
on both surface sides are terminated by hydrogen atoms. 
Geometry optimization is performed for the three top- 
most silicon layers until residual forces are less than 0.004 
Ry/A. The integration over the surface Brillouin zone is 
implemented using uniform 1 x8 k points. The images on 
the plane at about 3.1 A above the topmost silicon-atom 
layer are generated in the Tersoff-Hamann approximation 
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silicon-atom layer for several energy ranges, respectively. 
The characteristic peaks of the calculated LDOS agree 
with those obtained by the experiment Q. The filled- 
and empty-state images are mostly affected by the spa- 
tially localized bonding state and nonlocalized antibond- 
ing state, respectively, and the peak at 1.25 eV (1.0 
eV) below (above) the Fermi level Ep corresponds to 
the Si-Si and Si-H bonding (antibonding) states. Fig- 
ure ^ shows the simulated filled- and empty-state im- 
ages of various hydrogen-terminated Si(110) surfaces. In 
the filled-state images, the zigzag chains are uniformly 
oriented along the [110] direction, and neighboring hy- 
drogen atoms of two adjacent monohydride rows form 
a zigzag chain [Fig. Ela)]. The hydrogen atoms exist- 
ing at step edges form single-row chains of bright oval 
spots [Fig. EJL)], and an isolated row on a terrace is 
composed of a double-row chain of bright protrusions 
[Fig. 0£c)]. A missing row consisting of surface silicon 
and hydrogen atoms is observed to be a dark trough, 
and its edges are shaped by the single-row chains of 
bright spots [Fig. Efd)]. These images are in agreement 
with those obtained by the experiment £| . On the other 
hand, the empty-state images appear different from the 
filled-state images. The zigzag chains exist on the silicon- 
atom rows of the topmost layer [Figs. 0fe)-(h)] contrary 
to the filled-state images in which they are formed be- 
tween the two adjacent monohydride rows. The mono- 
hydrides existing at the step edges appear darker than 
those on the topmost layer [Fig.^Jh)]. The width of the 
dark trough composed of the missing row is observed to 
be two times greater than that in the filled-state image, 
and the contribution of the isolated row to the images is 
small [Fig. Hg)]. 

The difference between the filled- and empty-state im- 
ages is due to the contributions of the Si-Si bonding and 
antibonding states involved in tunneling: the filled-state 



Figures and [21 show the calculated LDOS of the 
atomically flat hydrogen-terminated Si(110) surface and 
the contour plot of charge distributions on the topmost 



FIG. 2: Calculated LDOS on the plane at 3.1 A above the 
topmost silicon-atom layer of the atomically flat hydrogen- 
terminated Si(llO) surface. The zero of energy is taken to 
the Fermi level Ef- 
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FIG. 3: (Color) Contour plots of charge distributions of LDOS for several energy ranges: (a) Ep — 1.5 to Ef — 2.0, (b) Ef to 
Ef — 1.5, (c) Ef to Ef + 1-1, and (d) Ef + 1-1 to Ef+1-5 eV. The contour maps are along the section in the (110) plane (left 
panel) and the (110) plane (right panel). The maps are illustrated on a logarithmic scale, and each contour represents twice 
or half the density of the adjacent contour lines. Pink and blue balls represent H and Si atoms, respectively, and large (small) 
balls indicate atomic positions on (below) the cross section. 



images primarily reflect the spatial distribution of the 
occupied bonding state, while the empty-state images 
are affected by the unoccupied antibonding state. In 
the fflled-state images, the positions between the silicon 
atoms of the topmost layer are dark troughs owing to the 
small contribution of the spatially localized Si-Si bond- 
ing state; therefore, the zigzag chains are observed to be 
formed by hydrogen atoms of two adjacent monohydride 
rows. On the other hand, since the Si-Si antibonding 
state is markedly expanded, the positions between the 
silicon atoms appear bright and the zigzag chains are 
formed on the surface silicon-atom rows. The monohy- 
drides existing at the step edges and in the isolated row 
appear darker than those in the terrace because the Si-H 
antibonding state of the monohydride at the step edges 
or in the isolated row shifts higher by ~0.6 eV due to the 
absence of neighboring monohydride rows. 

In summary, we have studied the electronic structures 
of various hydrogen-terminated Si(110) surfaces. The 
characteristic peaks of the calculated LDOS are con- 
sistent with the experimental ones. The formations of 
zigzag chains and bright protrusions in the filled-state 
images are in agreement with those obtained by the ex- 
periment On the other hand, in the empty-state im- 



ages, zigzag chains shift so as to locate on the silicon- 
atom rows of the topmost layer. The monohydrides ex- 
isting at the step edges appear darker than those in the 
terrace, and those in the isolated row hardly affect the 
empty-state images. The difference between these images 
is explained using the contributions of the bonding and 
antibonding states. To the best of our knowledge, this is 
the first study that indicates the differences between the 
filled- and empty-state STM images of the atomically flat 
hydrogen-terminated Si(110) surface, and an STM study 
is in progress to verify the results of this investigation. 
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FIG. 4: (Color) Simulated STM images of H-terminated Si(llO) surfaces: filled-state images of (a) terrace, (b) step, (c) isolated 
row, and (d) missing row models and empty-state images of (e) terrace, (f) step, (g) isolated row, and (h) missing row models. 
The filled-state images are obtained by integrating charge distributions from Ef — 2.0 to Ef eV, and the empty-state images, 
from Ef to Ef+1-5 eV. The meanings of the symbols are the same as those in Fig. The equi-LDOS surfaces at a height of 
about 3.1 A from the top-most silicon-atom layers are presented. 
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